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ABSTRACT 

The X-ray transient 080109, associated with SN 2008D, can be attributed 
to the shock breakout emission from a normal Type Ib/c supernova. If the ob- 
served emission is interpreted as thermal emission, the temperature and radiated 
energy are close to expectations, considering that scattering dominates absorp- 
tion processes so that spectrum formation occurs deep within the photosphere. 
The X-ray emission observed at ~ 10 days is attributed to inverse Compton 
scattering of photospheric photons with relativistic electrons produced in the in- 
teraction of the supernova with the progenitor wind. A simple model for the 
optical/ultraviolet emission from shock breakout is developed and applied to SN 
1987A, SN 1999ex, SN 2008D, and SN 2006aj, all of which have optical emis- 
sion observed at t ~ 1 day. The emission from the first three can plausibly be 
attributed to shock breakout emission. The photospheric temperature is most 
sensitive to the radius of the progenitor star core and the radii in these cases 
are in line with expectations from stellar evolution. The early optical/ultraviolet 
observations of SN 2006aj cannot be accommodated by a nonrelativistic shock 
breakout model in a straightforward way. 

Subject headings: shock waves — supernovae: general — supernovae: individual: 
SN 2008D 



1. INTRODUCTION 

The shock wave generated by the central explosion in a supernova accelerates through 
the outer layers of the star, giving rise to the peak photospheric temperature at the time of 
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shock breakout. The resulting X-ray /ultrav iolet flash was predicted long ago ( IKlein fe Chevalier 
19781 : iFalkl Il978t iMatzner fe McKee! Il999h. GRB 060218, discovered by its X-ray radi 



ation (|Campana et al. 2006 ). was associated with the fairly normal Type Ic SN 2006aj 



( Mazzali et al. 2006^. and its early emission has been interpreted as due to shock breakout 



( jCampana et al.l 120061 ; IWaxman et al.l 120071 ) . However, mildly relativistic motion is needed 



and it is likely that a central engine played a role in the emission. 

The X-ray transient XRT 080 109 was serendipitously discovered during Swift observa- 
tions of SN 2007uy in NGC 2770 (ISoderberg et al.l 120081 ) . The supernova associated with 
the burst, SN 2008D , was found to be spectroscopically similar to the Type lb SN 2005hg 
(IModjaz et al.ll2008l ). Early radio emission was detected, but not at the high luminosity that 
would suggest relativis tic motion, and no gamma ray burst emission was detected by Swift 
(ISoderberg et al.ll2008l ). In view of these properties, we further examine the possibility that 
the early (< 1 day ) emission from S N 200 8D w as entirely due t o the supernov a. This view 
was advocated by ISoderberg et al.l (120081 ). but Ku et al.l (120081 ) and |Lj ( 120081 ) find that a 
central engine flow is needed. We treat the X-ray emission in § 2 and the optical/ultraviolet 
emission in § 3. We take the distance to NGC 2770 to be 27 Mpc. 



2. EARLY X-RAY EMISSION 



The early X-ray emission from a supernova depends on the properties of the shock wave 
generated by the central explosion. The shock wave moving through the star is strongly ra- 
diation dominated, giving rise to a burst of radiation when the shock wave reaches an optical 
depth ~ c/u, where c is the speed of light and u is the shock velocity. The radiation accel- 
erates the gas ahead of the shock wave, so that the shock wave disappears; a gas dominated 
shock wa ve then forms in the surround ing medium when the supernova radius has roughly 
doubled ( Epsteinlll98ll ; iFranssonl 119821 ). This evolution is shown well in the simulation by 
Ensmanl (Il994f ). which has parameters intended for SN 1987A. The shock breakout begins 
at t m 5.8 x 10 3 s when the radiation shock starts to spread, and reaches a peak luminosity 
at t ~ 1.0 x 10 4 s when no shock is present. The gas shock starts to develop at t ~ 1.1 x 10 4 
s, when the luminosity is down by a factor of 3 from its peak value, and is well developed 
by t ~ 1.7 x 10 4 s, when the luminosity is down by another factor of 3. The simulation as- 
sumes pure absorption and an r~ 7 density power law for the preshock atmosphere, but these 
properties do not appear to be crucial for the shock evolution. Particle acceleration at the 
gas dominated shock gives rise to the radio synchr otron radiation emission that is observed 
from Type Ib/c supernovae, including SN 2008D. I Chevalier fe Franssonl (120061 ) found that 
for the conditions typically found in SNe Ib/c, the X-ray emission from the later interaction 
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is dominated by nonthermal mechanisms (inverse Compton and synchrotron) as opposed to 
thermal emission. 

Swift and Chandra observations of the X RT 080109 showed an e arly bright phase lasting 
~ 600 s, followed by a slow, fainter decline (jSoderberg et al.l 120081 ) . We identify the early 
bright phase with the breakout of the radiation dominated shock and the later emission 
with the interaction with the circumstellar medium. An important issue is the spectrum 
of the initial burst. There is general agreement that a power law spectru m (with photon 
i ndex T = 2.3) provide s a better fit than a blackbod y (with kT = 0.73 keV) (jXu et al.l 12008 



Soderberg et al 



2008 



Li 



2008 



acceptable (jXu et al. 2008 ; Modjaz et al. 



Modiaz et al.l 120081). although a b lackb ody fit may still be 



20081 ). ISoderberg et al.l (120081 ) find an absorption 



corrected fluence of 1.7 x 10 46 ergs for the power law spectral fit. 

Previous studies of shock breakout ( Klein Sz Chevalier 1978 ; Falkill978 ; Matzner Sz McKee 



19991 ) assumed that the initial breakout radiation has a blackbody or dilute blackbody spec- 
trum. The fits to the observed spectrum of SN 2008D show a preferenc e for a power law 



spectrum, and photon Fermi acceleration is invoked by Soderb erg et al. 



(}2008|) to explain 



the observations of SN 2008D, based on calculations by lWang et al.l (120071 ). In this scenario, 
the photons gain energy by scattering across a semi-relativistic shock front. The authors 
assume that an infinitesimally thin shock front is present, apparently from the time that 
the radiation dominated shock front breaks down at an optical depth of c/u. As discussed 
above, the breakout radiation is capable of accelerating the matter ahead of the shock front 
so that the formation of the gas dominated shock is delayed; the viscous shock starts to form 
only when the breakout luminosity is dropping from its peak value and the optical depth at 
the shock is decreasing. 

In view of the possible difficulty wit h producing a power law sp ectrum, we consider the 
thermal interpretation of the spectrum. iMatzner fc McKed (119991 ) estimated the breakout 
radiation expected from a relatively compact star with a radiative envelope; they determine 
the temperature by setting the postshock pressure = |aT 4 , where a is the radiation constant. 
For the explosion parameters, we use the energy E$i = 3, where E^i is the energ y in units 
of 10 51 ergs, and ejecta mass M = 4 M@, determined by ISoderberg et al.l ( 120081 ) from the 
supernova light curve and spectra. Taking the opacity, k, of the He rich gas to be 0.2 cm 2 
g _1 , the radius of the progenitor star required to produce kT = 0.73 keV at shock breakout 
is R — 1.5 x 10 11 cm. With these values of E, M, and R, the exp ecte d radiated energy at 
shock breakout is 1.7 x 10 45 ergs. However, iKlein fc Chevalier! (119781 ) and lEnsman fc Burrows 
( 119921 ) find that the observed color temperature can be 2 — 3 times the effective temperature 
if a blackbody were assumed because the spectrum is formed at a relatively deep layer in 
the scattering atmosphere. If the effective temperature is reduced by a factor of 2 to 0.36 
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keV, the radius is increased to R 
energy of 1.9 x 10 46 



ergs. 



6.4 x 10 11 c m = 9 R Q , which corresponds to a breakout 



Woosley et al.l (119951 ) find that a He star with initial mass 4 M ( 



© 



has a radius of 9 R Q at the time of carbon ignition; a He star can become extended in its 
late burning phases. With this radius, the peak shock velocity has ftj ~ 0.6, a high, but 
nonrelativistic, velocity. 

Although the thermal shock breakout model can reproduce the temperature and radiated 
energy of the observed burst, it falls short on the timescale. The expected timescale is 
~ R/c = 5 — 20 s, but the observed timescale is ~ 100 s (taken as the time with flux 
above half the peak flux). One possibil ity is that the photos phere is actually in a dense 
stellar wind around the Wolf-Rayet star (jSoderberg et al.ll2008l ). Assuming a constant wind 
velocity the optical depth to electron scattering in the wind is r e = 0.2/i e M_ 5 / (r 12 v 3 ) where 
the mean molecular weight fi e 0.5 for a non-hydrogenic gas, r 12 is the radius in units 
of 10 12 cm, the wind mass loss rate M_5 in 10 -5 M yr _1 and t>3 is the wind velocity in 
10 3 km s " 1 , normalized to a typical value for Wolf-Rayet stars. From the observed radio 



emission. 



Chevalier fc Franssonl (120061 ) estimated M-^v 3 1 = 2.8 for SN 19941 for an efficiency 



of magnetic field production in the shock front of 10%; a lower efficiency w ould increase the 



wind density. A similar estimate for SN 2008D, using radio data from ISoderberg et al. 



( 20081 ). gives M^v^ 1 = 1.5. For these parameters, we have r e = 1 at r 12 = 0.15, which is 
smaller than the estimated stellar radius, indicating that the freely expanding wind is not 
optically thick. 

However, the wind structure in the immediate vicinity of the star is likely to be more 
complicated than a p oc r~ 2 wind. The radiative acceleration thou ght to be responsible 
for th e h igh v elocity occurs over a region of seve ral stellar radii (e.g., 



Grafener &: Hamann 



20051 ; iLil 120071 ) , leading to an enhanced density. iGrafener fc Hamann! (120051 ) find that the 
acceleration up to ~ 1000 km s _1 takes place over a factor of 2 in radius, while the final 
velocity of ~ 2000 km s _1 was only reached at more than 10 stellar radii . Another important. 



but le ss understood, effect is connected to pulsational instabilities. iFadeyev &: Novikova 
(J200J) find that for He core masses below 10 M and L/M > 7 x 1O 3 L /M the star is 
pulsationally unstable. Hydrodynamical calculations show that these pulsations steepen to 
shock waves at the surface, which levitates the outer regions. For a 10 M He core and a 
luminosity of 2 x 1O 5 L they find that the radius of the star expands by a factor of two, 
and the scale height at the surface increases from 6 x 10~ 3 i? p h to ~ 0.23-R p h- There is 
observational evi dence for transien t mass loss before some Type Ib/c and Iln supernovae, 



e.g., SN 2006jc (IFoley et al.l 120071 ) . on a scale of ~ 10 16 cm. In view of these issues and 



the uncertainty in the nonequilibrium radiation field, the difference between the model and 
observed timescales may not be significant. 
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A crucial aspect of our scenario is that the spectral formation and thermalization occurs 
at a considerable depth in the supernova, so that the photons just scatter in the mass loss 
region. The place of thermalization is characterized by r scatt r abs « 1, where r scatt is the 
scattering optical depth and r abs is the absorption optical depth. We have r abs / r scatt « k s /k t , 
where the absorption is dominated by the free-free opacity Kg = 1.77 x l§~ 2 gn 2 T~ l ^v~ 2 in 
cgs units, kt is the Thompson opacity, and n is the density. The probability for absorption 
becomes r abs /r scatt « 1.46 x 10~ 23 gnT^ 1 ' 5 E^ / , where T 7 is the temperature in units of 10 7 K, 
.Ekev is the photon energy in keV, and cgs units are used. At 1 keV thermalization requires 
n > 6.8 x 10 22 /r s 2 catt cm -3 ; photons will only scatter at lower densities and thermalization 
occurs in the steep density region. We thus find that it is plausible that the temperature and 
energy can be reproduced, althoug h detailed hydro -dynamic , radi ative transfer calculations 
are needed to show this definitively. Ku et al.l ( 120081 ) andlLj (120081 ) argued against a thermal 
model for the emission based on the small radiating area derived from the luminosity and 
temperature. Here we argue that the spectrum is diluted in the scattering atmosphere, so 
that a larger radiating region can be accommodated. 



In a model intended for SN 19941, I Chevalier fc Franssonl (120061 ) found that inverse 
Compton emission dominates over days 1 — 20, with an X-ray luminosity on day 1 of 
vL v ~ 1 x 10 39 erg s _1 . The model luminosity evolves slowly because of the balance be- 
tween the shock wave expansion and the increasing optical luminosity, and on day 10 is 
~ 4 x 10 38 erg s _1 . The observed X-ray emission from SN 2008D by Ch andra on day 10 had 
a luminosity (0.5 — 10 keV) of ~ 10 39 erg s -1 (ISoderberg et al.l 120081 ). which corresponds 
to vL v ~ 3 x 10 38 erg s _1 for a L v oc v~ l spectrum. There is reasonable agreement with 
the model expectations; the radio and X-ray emi ssion suggest wind con ditions similar to 



those in SN 19941 on a scale ~ (2 — 3) x 10 15 cm. ISoderberg et ajj (120081) a lso advocate an 



inverse Compton origin for this emission. The I Chevalier fc Franssonl (120061 ) model predicts 
a drop in X-ray luminosity after day 10 when the photospheric lumi nosity drops. However , 
the supernova was close to the detection limit of Chandra on day 10 (ISoderberg et al.ll2008l ). 
making it difficult to check this prediction. 



3. EARLY OPTICAL EMISSION 



As the supernova photosphere expands and cools, the emission at lower frequencies, 
including optical, increases. Optical emission related to shock breakout thus peaks at a 
l ater time, ~ a day or day s. Such emission was seen in the early observations of SN 1987A 



(IHamuy et al.l Il988l ) and IChevalierl (119921 ) developed an analytical theory to describe the 
emission. The main approximation in this approach was that the supernova be well into its 
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free expansion phase, which requires a relatively small radius progenitor star. The model 
should thus be applicab le to Type Ib/c s upernovae and we consider the Type Ic's S N 1999ex 
( Istritzinger et al.l l2002h and SN 2006aj JCampana et alibood : iModiaz et aDbood ). as well 
as SN 2008D. In a ll these cases, there are observations at an age ~ 1 day. Here, the theory 
of I Chevalier! ( 119921 ) is extended so it can be compared to detailed observations. 



After correcting for a typographical error, equation (2.8) in IChevalierl ( 119921 ) for the 
density of outer freely expanding gas at velocity v becomes 



p = 1.4 x 10 98 £fi 59 M 



2 "(I/,. 7.18^-3^-10.18 g ( . n| 3 



) F 7. 



X 



where M is the ejecta mass in units of M & and F\ is a factor that describes by how much 
the velocity of an element of gas increases after it has first been shocked. This expression 
was derived for the case of SN 1987A, which had a progenitor with a radiative envelope, like 
SNe Ib/c. In a planar blow off F\ = 1.81, but is reduced in a spherical expansion. Taking 
■ Fi = 1.35 leads to density and pressure profiles that agree with results of iMatzner &: McKee 



(119991 . their equations 46 and 47) for the outer parts of an exploded star with a radiative 
envelope, and we use that clS 8u reference value. 

Given the density structure, the photospheric radius can be found by integrating np in 
to the point where r = 1. Taking k = 0.2 g cm -2 , the photospheric radius is then 



r p h 



3.0 x 10 i4 E 5 V 9 M - a ^(F 1 /1.35) 



°' 78 # 78 cm, 



where td is the age in days. The velocity of freely expanding gas at the photosphere is 



Vph 



3.4 x lO^fM- ' 28 



(Fi/1.35) - 78 * 



a7 V 22 km s" 1 . 



(2) 



(3) 



The luminosity at the phot osphere can be f ound by allowing for a diffusion wave to move 
back into the exploded star (|Chevalierill992l ) 



3.3 x 10 42 £ 5 °i 91 M -°- 74 



/^(Fi/l.SS)- 017 ^ - 34 erg s" 1 , 



(4) 



where R\2 is the radius of the progenitor star in units of 10 12 cm. Because the inner layers are 
shocked at a lower velocity and are observed at a later time than the X-ray observations, the 
radiation field in this case can be approximated by a blackbody. Using L c = 47rrp h aT^, where 
cr is Stefan-Boltzmann's constant and T e is the effective temperature of the photosphere, we 
have from equations (j5J) and (j3J 



1.4 x 10^ 5 Y 3 M -° U4 i^f (Fi/1.35) 



-0.43^-0.48 K 



(5) 



These results can be compared with those of IWaxman et al.l (120071 ). who also assumed 
that the photosphere is in the outer shock accelerated part of the supernova density profile, 



-7- 



but otherwise used a different method. The scaling with para meters is similar , altho ugh not 
exactly that found here, and, for the same reference values, I Waxman et al.l (120071 ) obtain 
r ph — 3.2 x 10 14 cm and T e = 2.6 x 10 4 . The agreement with the r p h found here is good, but 
the value of T e found here is smaller by a factor of 1.9, or the luminosity is smaller by 12. 
The method of Waxman et al.l (120071 ) does not include radiative diffusion, whereas that of 
Chevalierl (119921 ). used here, does. We find that the luminosity decreases by 1.9 because of 
diffusion, for the standard parameters at = 1, a relatively small effect. 

The photospheric radius and temperature determine the spectral luminosity evolution, 
as shown in Figure 1 for 4 wavelengths relevant to the Swift UVOT bands. The results in 
Figure 1 assume the reference values for the physical parameters (E 5 i — 1, M — 1 M , i? 12 = 



1), but can be used for other values by scaling the luminosity by a factor E®f 8 M Q 0S8 R®f 1 



and the time by a factor E^^ M^R^ . When the wavelength is in the Ray leigh- Jeans 
part of the spectrum, we have L v oc t 



i.i 



For the theory described here to be applicable, the photosphere must be within the 
steep power part of the supernova density profile. If there is no surrounding medium, the 
limitation at the high velocity end is th e maximum velocity at s hock breakout, which is 



129, OOOE^Mo" - 42 ^ - 32 km s" 1 (IMatzner fc McKeelll999h . Using equation © and 



the reference values, this velocity corresponds to the velocity at the photosphere at an age of 
0.002 day. If the supernova is surrounded by a dense, optically thick wind, the time at which 
the theory becomes applicable is after the shock wave has broken out from the wind. Also, 
because the density distribution at small radii becomes flatter than the steep power law, 
the actual value of r p h should gradually become smaller than the value estimated here. To 
estimate this effect, we used the harmonic mean density profile of IMatzner fc McKed (119991 ) 
(their eq. [46] ) to estimate the value of r p h more accurately. For the reference parameters and 
an age of 1 day, the value of r p h is 20% smaller with the more accurate density distribution. 

Another requirement of the applicability of the simple model is that the opacity remain 
constant. Recombination of He ++ to He + occurs at a temperature of a few 10 4 K at the 
densities of interest; this results in a decrease in the opacity by a factor of 2. A larger 
decrease in opacity occurs when the gas becomes neutral, which occurs at T ~ 1.6 x 10 4 K. 
The recombination adds internal energy to the gas, which can help power the light curve. 
These considerations show that the simple model for the early light curve starts to break 
down at an age ~ 1 day. 



The earliest observations of SN 1987A by lHamuy et al.l (119881 ) are on day 1.14 when 
L = 4.9 x 10 41 erg s _1 , r ph = 1.4 x 10 14 cm, and T e = 13, 600 K. To apply the model developed 
here, we use the parameters k = 0.34 cm 2 g _1 (to account for a H and He composition), 



Ei 



51 



1, M = 16, and R 12 = 2.1 with the result L = 5.5 x 10 41 erg s" 



r P h 



1.6 x 10 



ii 
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cm, and T e = 1 3, 100 K. In this case, the par ameters are fairly well determined by extensive 



modeling (e.g., Ilmshennik &: Nadezhinlll989l ) and it appears that the model presented here 



can approximately represent the early evolution. 

There are good early limits on the brightness of the Type Ic SN 1999ex because of the 
occurrence of another supernova, SN 1999ee, in the galaxy, implying that the age of the 
supernova was < 1.05 days when first obse rved. At that time, L = 2.8 x 10 41 erg s~\ r ph = 
1.2 x 10 14 cm, and T e = 12, 900 ± 1400 K Jstritzinger et al.ll2002h Istritzinger et all f|2002h 
compare the light curve of SN 1999ex to models of IWoosley et al.l (119871 ) with M = 6.2 M ; 
with E51 = 1, our model then gives r p h ~ 1.8 x 10 14 t°' 78 cm, which is in adequate agreement 
with the observations. The corresponding model value of T e is 13, OOOi?^ 25 ^ 8 K. 



From Swift UVOT observations of SN 2008D, ISoderberg et al.l (120081 ) estimate that 
r p h ^ 3 x 10 14 cm and T w 10 4 K at an age of ~ 1 day; the closest time of observation 
is 1.7 day, so these values must be regarded as estimates. For E51 = 3, M = 4 M & , and 
R = 9 Rq, our model gives r p h = 3.1 x 10 14 and T = 1.2 x 10 4 , in reasonable agreement with 
the observed values. We note that the radius R refers to the outer star layers where shock 
acceleration and a steep density gradient occur. There is the expectation in the model that 
the supernov a should be cooling at this age, which is consistent with the optical/ultraviolet 
observations ( Soderberg et al.l I2OO8I ) . 

For SN 2006aj, ICampana et al. I (bood ) estimate r p h « 3 x 10 cm and a temperature 



of 3 eV at t = 10 5 s = 1.16 day, corresponding to T e w 35, 000 K and L fh 1 x 10 44 erg s _1 , 
far more luminous than either SN 1987A or SN 1999ex at a comparable age. To compare 
the observations with models, we use E51 = 2 and Mq = 2 (IMazzali et al.ll2006l ). The model 
reference properties on day 1.16 values are then r p h = 3.6 x 10 14 cm and T e = 13,000 K, 
corresponding to L = 2.6 x 10 42 erg s _1 . The value of r p h is in reasonable agreement with 
observations, but the model T e is low, leading to a significantly lower luminosity than is 
observed. The critical parameter for T e is the radius of the progenitor star. The required 
radius to produce the temperature is R ~ 5 x 10 13 cm, which is the radius of an extended 
red supergiant star and much larger than that expected for a Wolf-Rayet star. As noted 
above, a helium star can become more exten ded in the evolution leading up to the explosion 
if its initial mass is low (IWoosley et al.lll995l ). but the effect does not appear to be sufficient 
to explain the radius needed here. 

In addition, the evolutio n leading up to 10 5 s can be considered. The Swift UVOT 
light curves for GRB 060218 (ICampana et al.l I20Q6L see also Ghisellini et al. 2007) can be 
compared to the model curves in Fig. 1. The cooling indicated by the evolution of the 1880 
and 2510 A model light curves matches that seen in the observations fairly well. However, 
between 10 4 and 10 5 s, the observations show a turnover at the longer wavelengths (3450 and 
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5440 A) that is not in the model curves and cannot be attributed to the expected cooling of 
the photosphere. The turnover indicates a decrease in the emitting area that is not expected 
if the emission is from the outer part of the exploding star. In addition, it can be seen that, 
although there is a roughly power law increase in the fluxes at early times, the observed 
slope is shallower than the t 11 dependence expected in the model. For these reasons, the 
prob lem of underproductio n of flux in the model at 10 5 s is even more severe at early times 



sec 



Ghisellini et all 120071 ). Although the early Swift UVOT observations of GRB 060218 
show a tantalizing similarity to expectations for early supernova emission, a shock breakout 
model does not reproduce the observations in a straightforward way, indicating that a central 
engine plays a role in the emission. 
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Fig. 1. — Model predictions for the photospheric spectral luminosity in the shock breakout 
phase. The curves are for 1880 A (dash), 2510 A (dot), 3450 A (long-dash), and 5440 A 
(dash-dot). The model assumes reference values of the physical parameters, but can be 
scaled to other values, as described in the text. 



